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ABSTRACT 

In order to be in a long-lived configuration, the density in a fluid disk should be constant 
along streamlines to prevent compressional (PdV) work from being done cyclically around every 
orbit. In a pure Kepler potential, flow along aligned, elliptical streamlines of constant eccentricity 
will satisfy this condition. For most density profiles, differential precession driven by the pressure 
gradient will destroy the alignment; however, in the razor-thin approximation there is a family 
of simple equilibria in which the precession frequency is the same at all radii. These disks may 
therefore be long-lived at significant eccentricities. The density can be made axisymmetric as 
r —^ 0, while maintaining the precession rate, by relaxing the requirement of constancy along 
streamlines in an arbitrarily small transition region near the center. In the limit of small eccen- 
tricity, the models can be seen as acoustically perturbed axisymmetric disks, and the precession 
rate is shown to agree with linear theory. The perturbation is a traveling wave similar to an 
ocean wave, with the fluid rising and falling epicyclically in the gravitational field of the central 
mass. The expected emission line profiles from the eccentric disks are shown to be strongly 
asymmetric in general, and, in extreme cases, prone to misinterpretation as single narrow lines 
with significant velocity offsets. 

Subject headings: galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: structure — 
stars: formation — stars: novae, cataclysmic variables 



1. Introduction 



Accretion disks, protostellar disks, and galactic disks are almost always assumed to be axisymmetric 
or nearly so. This assumption is not without justification. Disks support a rich spectrum of stable and 
unstable modes which can efficiently redistribute angular momentum while viscous, collisional, and radiative 
processes dissipate energy. Circularization is then a naturally expected result. Nonetheless, disks of finite 
eccentricity may play an important role in certain astrophysical systems, and, indeed, appear to be required 
by observed phenomena on scales from AU to parsecs. 

Two examples illustrate this need particularly well. At the stellar scale, eccentric disks are well estab- 
lished as the source of the superhump phenomenon in short-period cataclysmic variables. The period of these 
transient luminosity variations differs by a few percent from the binary orbital period, and the difference 
is readily explained by a precession of the disk (Vogt 1982; Osaki 1985). The eccentricity is "pumped" 
by the binary tidal field; numerical and analytic studies (Whitehurst 1988; Lubow 1991) indicate that the 
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eccentricity can be excited at the 3:1 Lindblad resonance (but see Heemskerk (1994) and Stehle (1999) for 
dissenting views). While the eccentricity may vary between outbursts, the repeatabihty of superhump light 
curves within an outburst (Patterson 1998) indicates that the eccentricity is not washed away on short time 
scales by differential precession. Instead, the disk needs to precess coherently for at least several orbital 
periods. The primary driver for this precession is most likely the gravity of the secondary star, but the disk's 
own pressure gradient may also bo a significant influence (Murray 2000). 

At the parsec scale, the best studied eccentric disk candidate is the "double nucleus" of M31 (Lauer 
et al. 1993; Tremaine 1995; Kormendy & Bender 1999). This disk is stellar rather than gaseous, and is 
essentially isolated in the Keplerian potential of the central black hole, since neither the bulge nor the main 
disk has significant influence at this scale. Assuming the disk is not a transient, its structure must be such 
that its own self gravity is able to drive a coherent precession. This condition imposes a characteristic non- 
monotonic radial eccentricity profile on the disk (Statler 1999; Salow & Statler 2001), which is reproduced 
by dynamical simulations (Bacon et al. 2001). However, whether the eccentricity can be self-excited or must 
be tidally driven remains at issue (Bacon et al. 2001; Tremaine 2001), and the dynamical stability of such 
conflgurations is midctermined. 

This paper takes up a question complementary to both of the above examples, namely, whether iso- 
lated Keplerian fluid disks can exist in long-lived eccentric equilibria. I will not address directly issues of 
hydrodynamic stability, which are extremely complex — even in axisymmetric systems — and far beyond the 

simple arguments presented here. But some necessary constraints on disk structure can be gleaned without 
a lengthy stability analysis. A slowly evolving disk should have a density that is approximately constant 
around streamlines. Were this not the case, PdV work would be done cyclically, and most likely irreversibly, 
aromid every orbit, dissipating the eccentric motions. In addition, the density and pressure must be such 
that perturbed streamlines can precess coherently. In the approximation of a razor-thin two dimensional 
disk with a polytropic equation of state and negligible mass, these minimal criteria for astrophysical realism 
are sufficient to define a remarkably simple family of eccentric disk models. 

The arguments of the paper are laid out as follows: Section 2 presents the basic properties of the models. 
Simple scaling arguments are used in § 2.1 to derive the radial density profile for an assumed equation of 
state. Section 2.2 calculates the pressure-driven precession, and demonstrates that in the frame that rotates 
at the precession frequency, the motion of fiuid elements is Keplerian to first order in the pressure. In § 
2.3 the slightly sticky question of the inner boundary condition is considered. It is shown how the central 
eccentricity can be taken smoothly to zero by altering the density profile in an arbitrarily small transition 
region. The structure of the models in the low-eccentricity limit is shown in § 3 to correspond to that 
of axisymmetric disks perturbed by an m = 1 traveling wave. The expected emission-line profiles from 
unresolved disks are computed in § 4. Finally, § 5 discusses the connection with other work on eccentric 
disks, describes the possible instabilities that could affect the models, and reiterates the main results. 

2. Disk Models 

2.1. Scaling Arguments and Radial Profile 

The disks of this paper are assumed to be planar and razor-thin, with zero thickness perpendicular to 
the plane. The fiuid has zero viscosity, and is characterized by a surface density S and a two-dimensional 
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pressure 11, related by a 2-D analogue of a polytropic equation of state, 

n = iiTE^+i/", (1) 

in which n is the polytropic index. Sound waves travel at a speed 

" 'n+ 1\ n 



c = 



1/2 

(2) 



and a fluid clement feels a forcc^ — VII/E in the presence of a pressure gradient. The force can also be 
written as the negative gradient of the enthalpy H, where 

if = ii:(n + l)Ei/". (3) 

The gravity of the disk is ignored. 

A disk of non-interacting particles (or pressure-free fluid elements) on aligned Kepler orbits has a surface 
density given by (Statler 1999) 

In equation (4), a is the scniiniajor axis and E is the eccentric anomaly. The arbitrary function /i(a) gives 
the mass per unit interval of a, and the eccentricity profile is described by e(a) and its first derivative e'{a). 
Equation (4) is valid only if the orbits do not cross, which requires |e + ae'\ < 1; but e is not assumed to be 
small. (A more general version of this formula for unaligned orbits is given in Appendix A.) 

Clearly if e' = 0, then S = 5](a) and the density is constant around each orbit. This happens because 
the speed at each point around the orbit is inversely proportional to the local separation between neigh- 
boring orbits, and so competing terms in the equation of continuity cancel. Now, treating the pressure as 
a perturbation, the equation of state (1) implies that the pressure is stratified on the unperturbed orbits. 
The perturbing force on each fluid element is therefore perpendicular to the orbit, and does no work. Its 
primary effect is to drive a precession, at a rate Ctp which must be constant at all radii for the disk to be 
long-lived at finite eccentricity. This requirement is sufficient to determine the radial density profile of the 
disk, as follows: The small dimensionless parameter in the problem is e = f2p/w, where oj is the Keplerian 
mean motion. Since the latter scales as a~^/^, ftp = constant implies e ~ a^/^. To leading order, e must 
be proportional to the ratio of the perturbing force — VII/E to the Keplerian force GM,/?-^, where M,. 
is the central mass. These quantities should each be thought of as averaged over the unperturbed orbit; 
since e = constant, these averages are proportional to —{dIl/da)/J^ and GMc/a?, respectively. The required 
scaling for e then implies that {dR/da)/11 ^ a~^/^, which can be written as 

;,(l + l)s-«/.,s = -(±)"'^.fc. (5) 

where a* is an arbitrary constant. Integrating this equation and defining ai = [K{n + l)/2{n — l)]^a» gives 
the density profile 

E(a) = So[l-(a/ai)i/2l"^ 



1 Throughout this paper "force" refers to force per unit mass. 
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where the surface density is equal to Sq at the center and falls to zero at a = ai . 

In the limit n ^ oo the equation of state (1) becomes isothermal and one obtains the density law 

S = Socxp [-(a//i)i/2l ^ 



(7) 



where /i is a constant scale length. However, this disk is infinite in extent, requiring e to become arbitrarily 
large since it must scale as c?!"^ . I therefore exclude the isothermal case as unphysical. 



2.2. Precession Rate and Surface Density 

Fluid elements in the disk move along perturbed Kepler orbits. Ordinarily one would expect, in addition 

to the secular precession, cyclic variations in the orbital elements that would distort the pure ellipses and 
necessitate a first order correction to the surface density. In this case, however, one can show that the 
correction is zero. 

Using equations (2), (3), and (6), the enthalpy can be written as 

H = ncl \l - [a/aif/^] , (8) 

where co is the sound speed at a = 0. The radial and tangential components, R and T, of the perturbing 
force are therefore given by 

T^ _ f ^ 1^ 

^ ' ' da \dr'rdf 

2 

where / is the true anomaly and a = r(l + ecos/)/(l — e^). But the Keplerian velocity has components 

= ^Y72(Y3^2yT72(esin/,l + ecos/); (10) 

consequently, in the frame that rotates at speed fip, the Coriolis term —2^pZ x v will identically cancel the 
pressure gradient if 

^ - ill) 

In this case, the disk is globally in geostrophic balance, and the only force term surviving in the Eiiler 
equation is the Kepler term —GM^/r^. (The centrifugal term n^r is second order.) Thus in the rotating 
frame the fluid elements follow exact Kepler ellipses, and the surface density is given by equation (6) with 
no first-order correction. 

It is easy to see intuitively why this result holds. Since the density is constant around each orbit, the 
perturbing force is proportional to the pressure gradient, which is inversely proportional to the separation 
between isobars. But the isobars are the orbits, whose separation is inversely proportional to the orbital 
speed when e = constant. Thus the perturbing force is both proportional to the velocity in magnitude and 
perpendicular to it in direction, just as the Coriolis force is. The two forces can cancel on all orbits for the 
same Op owing to the a^/^ dependence of the enthalpy. 
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Fig. 1. — (a) Surface density contours for the disk with n = 3 and e = 0.5. Contours are logarithmically spaced by factors 
of 2. (b) Surface density as a function of scmimajor axis for models with values of the polytropic index n as indicated, (c) 
Semimajor axis enclosing half the mass as a function of n. 



One can verify that equation (11) agrees with the precession rate obtained from standard perturbation 
theory by using equation (9) in the formula for the instantaneous drift rate of the argument of pericenter, 



1 

zu = — 

e 



'a{l-e^y 









-^cos/ + fsin/^±l^^) (12) 
^ 1 + ecos// 



(Burns 1976; Murray & Dermott 1999), and averaging around the orbit. Equation (11) may also be written 
in the form 



, 2 

n / Co 



"^ = -4(l-e2)i/2 UJ -1' (13) 

where vi and ui are the Keplerian speed and mean motion at the disk edge. This shows that the expansion 
parameter e is basically a''/^ times the square of the ratio of the central sound speed to the edge Keplerian 
speed. Notice also that the precession is retrograde, and is faster for larger eccentricities. 

Figure la shows a density contour plot for a disk with n = 3 and e = 0.5, illustrating the stratification 
of the density on orbits. Figure lb compares the radial density profiles for selected values of the index n. 
The total mass of the disk is 

M = 24.^1^(1 -e^)V^Eoa?, (14) 

where r(a;) denotes the gamma function. The semimajor axis ah enclosing half the mass is given by the 
solution of 

4(4,n + l) = ^, (15) 

where Ix{p, o) is the incomplete beta function and z = (a/j/ai)^/^. Figure Ic shows the variation of ah with 
n. For n = the density is constant and ah = 2~^/^ai; as n increases the disks become increasingly centrally 
concentrated. 



2.3. Behavior at the Center 



The behavior of the models near their centers presents a problem. It is not realistic to expect e to 
remain finite as o ^ 0, since there are no strongly eccentric central masses in astrophysical systems. One 
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could imagine carving out an inner hole from the disk to avoid the center; but this expedient would remove 
the pressure gradient needed to synchronize the precession of the innermost material with that of the rest of 
the disk. Instead, the eccentricity needs to go smoothly to zero at the center — or at small finite radii — while 
maintaining the precession rate. I show here how this transition may be effected by allowing the density 
to deviate from that in equation (6) over a small interval of a, in the approximation where the density and 
eccentricity changes are both small (in a sense defined below). Formally this assumption limits the treatment 
to eccentricities < 0.3; but it seems reasonable that a similar approach can be taken for larger e. 

Readers more interested in the results than in the technical details are invited to skip down to the last 
paragraph of this section. 

At leading order, the changes to the precession frequency from the altered density and eccentricity 

profiles can be considered separately and summed. I begin with the former, replacing the surface density S 
by E[l + 6{a)], where S is the fractional density increase. The equation of state gives 



vn 



(16) 



where, from equations (3) and (8), K = [n/{n + 1)]cqEq The first term on the right-hand side is 
responsible for the precession given in equation (11). Because 6 is constant on orbits, the second term 
simply contributes an additional d/n times the same rate. The third term is more difficult. The radial and 
tangential components of V5 are 



86 ^, 1 + e cos / 



dr 1 — r df 

where S' = dS/da. Using equations (6), (16), and the relations 



1 dS ^, esin/ 



1 



cos/ 



cos E — e 
1 — e cos E ' 



sin/ 



(1 - e^/^ sinE 
1 — e cos E ' 



this results in the force components 

{R,f)=cl\l-{a/a^)'/^ 



6' 



sinE 



(17) 



(18) 



(19) 



1 - ecosE' (1 - e2)V2(i _ ecosE) 

The instantaneous precession rate can be obtained from equation (12), written in terms of the eccentric 
anomaly using equation (18). Because udt = (1 — e cos E)dE, the resulting expression can be timc-avcragcd 
around the orbit by multiplying by (1 — ecosE) and averaging over E. (Useful formulae for evaluating the 
necessary averages are given in Appendix B.) Adding the contributions of both the second and third terms 
in equation (16), the total extra precession induced by the density enhancement 5 is found to be 



'p5 



2[GMcai(l-e2)]i/2 

js S'[{aa^)y'-a] 
I 2 e 



1 + 2 



(l_e2)(l-e) (l-[l-e2]i/2) 



(20) 



Note that Qps < for positive density enhancements that diminish outward. 



To calculate the corresponding result for the eccentricity gradient, S is replaced by S[l + g{r, /)], where, 
from equation (4), 



, 2e+ (l + e2)cos/ 
9ir,f)=er (j^-^ip • 



(21) 
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The equation of state gives, to leading order in g, 



(22) 



As in equation (16), the first term on the right hand side gives the ordinary precession rate (11). However, 
the second term does not produce a proportional contribution because g is not constant around the orbit. 
Instead, it produces the force components 



1 



cos_E + e esin£'(cos£' + e) 

1 - ecos£;' " (1 - e2)i/2(i _ ecos£;) 



The third term is algebraically cumbersome because of the gradient of g{r, /). Defining 



gi{a) = e{e' + ae") + 2a 



1 



93{a) 
9i{o) 

the components of Vg can be written as 

dg gi{a} + g2{a) cos E 



l-e2 
4aee'^ 



1-e 



2 ' 



egi{a) + e', 
6^2 (a) - e']. 



1 dg [93{0') + 54(0) COS E] sin E 



dr (l-e2)(l-ecos£;)' r df {1 - e'^f/'^ {1 - e cos E) ' 
This produces the force components 

e' r gi{a) + g2{a) cos E [gs{a) + g4{a) cos e] sin E 



{R,T) 



1 - (a/ai)i/2 



1 -e2 



1-ecosi; ' (l-e2)V2(i_ecosi;) 



(23) 



(24a) 

(24b) 

(24c) 
(24d) 



(25) 



(26) 



Following the same procedure as above, I find the total extra precession induced by the eccentricity gradient 
to be 



Co 



1 



(GMeai)V2 e(l - e^)^/^ 



ae 



+ 



(aai) 



1/2 



■^Q(a)gj(a) 



where the Cj coefficients are given by 

ci(a) 
C2(a) 
C3(a) 
C4(a) 



-e(l-e2) [l-(l-e2)i/2' 
(1 - e2)3/2 1 - (1 - e2)V2 

l-^-(l-e^)^^ 



2 ^ 



(27) 

(28a) 
(28b) 

(28c) 

(28d) 



The inner transition region, where e is varying and 6 is nonzero, should join smoothly onto the main 
body of the disk, where e = eo = constant. The requirement that both regions precess together is given by 



rip -j- -j- Opg/ — ? 



(29) 
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Fig. 2. — (a) Eccentricity e {thin lines) and fractional density enhancement S {thick lines) in the small inner region where e 
goes to zero according to equation (31), for models with n = 3. Successive curves, bottom to top, show results for main-body 
ellipticities eo = 0.1, 0.2, 0.3, 0.4, 0.5. (6) Surface density contours {solid lines) in the inner region of the model with n = 3 and 
eo = 0.3. Contours are logajrithmically spaced by factors of 2^/®. Dotted lines show representative fluid streamlines. 



where ftpo is the main body precession rate, i.e., equation (11) evaluated at e = eo- This requirement yields 
a differential equation linking e(a) and d{a): 



^[(aat)V2 - a] {e^ + 2(1 - e^)(l - e)[l - (1 - e^)V2]} 



= S + n 



1-e" 
1 - p2 

1 Cq 



- 1 



+ 



e(l - e2) 



-ae' + 4- 



4: 

-^c^(a).gj(a) 



(30) 



Since the point of the exercise is to manipulate e, I let e(a) be a fixed function and use equation (30) to 
calculate 6{a). In order to avoid a discontinuity in the gradient of the surface density, e(a), e'{a), and e"{a) 
should all be continuous at a = ao- A function with the desired properties is 



e(a) 




a < ao; 

ao < a < ai. 



(31) 



For < 6 < 1, e(a) is linear at small a and joins smoothly onto e = eo at a = ao- 



Figure 2a shows profiles of e and the fractional density enhancement S in a transition region with radius 
ao = O.Olai and b = 1/2, for a set of disks with n = 3 and various eccentricities Cq. The results are fairly 
insensitive to the value of n. One sees that the eccentricity can be made to go continuously to zero at the 
center by steepening the density profile through most of the transition region. Figure 2b shows density 
contours in the transition region for the disk with eq = 0.3. While the density is no longer constant along 
streamlines, the variation about the mean is in the worst case only ±20%. For eo > 0.3, the needed density 
enhancement becomes of order unity, which violates the assumption behind equation (16). The assumption 
that g{r, /) is also small requires that ae' ^ 1 — e, a condition that is consistent with, but more stringent 
than, that for the non-crossing of orbits. The choice of e(a) in equation (31) begins to violate this condition 
foi' Co ^ 0.3. But these conditions are required only by the first-order expansion; a more careful treatment 
may be able to find similar remedies for larger eo- 
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3. Connection with Axisymmetric Disks 



For small e, the structure of the models can also be understood in terms of the acoustic modes of an 
axisymmetric disk. A formalism for dealing with such disturbances is given by Papaloizou & Savonije (1991) 
and Hecmskerk, Papaloizou, & Savonije (1992, hereafter HPS). They start with a perturbation to the surface 
density of the form 

E = E(r) + j:'(^r)e'^rn4>+<rt) ^ ^32) 

where S(r) refers to the unperturbed disk; similar forms are assumed for the other hydrodynamic quantities. 
They then merge the linearized fluid equations into a single operator equation (HPS's eqs. [9] and [10]): 



S' = jC{W), 



(33a) 



where 



dW 2mnaW 
dr K^r 



'dW 2rnnaW 
dr K^r 

2n d{Pn) 

r dr ' 
1 d{r^n) 



2mil(TE niW d fl 



K'^rD 



ar dr \( 



a 
D 



rS dr 
a + mfl, 



_ -2 



(33b) 
(33c) 

(33d) 
(33e) 
(33f) 



and f2 is the circular frequency. In the simple case where the only gravity is that of the central mass, 
W = c^E'/E, and equation (33) (rather than HPS's more complicated equation [15]) governs the modes. 



In the axisymmetric limit, the present models have a surface density given by 

E(r) = Eo [l-(r/ai)V2]" 

and a sound speed given by 



(r-M) 



1/2 



The pressure gradient reduces fl from its Keplerian value to 



1/2 



A{GMcai 



u/2- 



(34) 
(35) 

(36) 



The difference in surface density between the eccentric and axisymmetric disks is, in the limit of small e. 



AE = Eo { [l - (r/ai)^/2 (l - |cos/)]" - [l - (r/ai)^/^]"} . 



(37) 



For arbitrary n, AE is a superposition of Fourier harmonics. I restrict the discussion here to the case n = 1, 
so that 



AE 



EqC 



1/2 



— (-) COS/. 



(38) 
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This corresponds to a single m = 1 term, with 



^' = ^[-]'\ (39) 



(40) 



and consequently 

W ='■ 

2 Voi, 

Direct substitution into equation (33b) shows, after lengthy algebra, that equation (33a) is satisfied to leading 
order in Cq when a = CQ/[4(G'Mcai)^/^]. The form of the exponential in equation (32) indicates that this 
value corresponds to a regression of the pattern and agrees with equation (11) for n = 1 and e = 0. Also, 
since a is real, the mode is stable. 

The reader may have noticed that the inner and outer radial boundary conditions have not played a role 
in this discussion. The fact that S' is real means that the mode has zero radial wavenumber. Unlike WKB 
waves, this wave is completely unwound and propagates only in the tangential direction. In this respect it is 
not really a "mode" in the sense of being a standing wave in a resonant cavity. Instead, it is a traveling wave. 
In fact, since the wave does not compress the fluid — even though the fluid itself is not incompressible — it 
bears a close similarity to an ocean wave. Like an ocean wave, fluid elements move epicyclically with respect 
to the mean flow, and the restoring force for the pressure driven vertical displacement is the gravity of the 
central mass. The variation of the sound speed and flow velocity with depth conspire to refract the wave 
around the central mass while keeping the wave front straight. 



4. Line Profiles 

Constant-eccentricity Keplerian disks are kinematically simple. Spatially resolved rotation curves along 
any line through the central mass will show velocities proportional to r^^/^, but will not, in general, be 
antisymmetric about the center. The observed velocities on opposite sides of the central mass will differ by a 
multiplicative factor depending on both the eccentricity and the viewing geometry. Kinematic measurements 
along several position angles would be necessary to determine the disk structure and orientation; but a 
reasonably accmatc estimate of the central mass could be obtained by averaging the masses derived from 
each side of the rotation curve under the naive assumption of circular motion. 

It seems more likely, however, that disks of the sort discussed here will usually be unresolved. The 
expected emission line proflles from an unresolved disk can be straightforwardly calculated. For simplicity, 
I assume the emission comes from a recombination line of a species whose density follows the total density, 
so that the emissivity is proportional to S^. The disk is assumed to be optically thin with the line of sight 
in the disk plane, an angle 9 from the major axis. A mass dM on a single Kepler orbit in the disk produces 
a contribution to the line of sight velocity distribution (LOSVD) f{vf ) given by 

... . ^ dM {1-ecosEf 

^^""^^ 27r(GM,/a)i/2 (cos i;-e) cos + (l-e2)i/2 sin E sine ■ ^ ' 

At each V£, the contribution must be summed over the two values of E given by 

. ^ A -5(^2 +52 -1)1/2 B + A(A2 +52-1)1/2 

-^^^1 = A2 + 52 ^osE, = ^2+^2 (42) 
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Fig. 3. Simulated line profiles for unresolved elliptic disks. Relative fiuxes per unit interval of velocity are shown for velocities 
in units of (GMc/ai)^/'^. Results for eccentricities of 0.1, 0.4, and 0.7 are shown in the first, second and third columns, and for 
lines of sight along, oblique to, and perpendicular to the disk major axis in the top, middle, and bottom rows. Different line 
styles correspond to values of the polytropic index n as indicated at top left. The central sound speed is 0.2. 



and 



where 



sini?2 = 



A + B{A^ + B^ - 1) 



1/2 



742 + 52 



cos £^2 = 



B - A{A^ + 52 _ 1) 

A2 + B2 



1/2 



A 
B 



{GMJaf/^ 



{GM,/aY'^ 



cos 6, 



(l-e2)V2sine + e. 



(43) 

(44a) 
(44b) 



These expressions ignore the precession of the disk, which adds a slow, retrograde, soHd-body component 

to the rotation; this effect is included in the examples shown below, but makes only a tiny difference to the 
results. Assuming an ideal gas, the equation of state and the density profile imply a one-dimensional thermal 
velocity dispersion 



2 



1 



1 - (a/ai)i/2 



(45) 



The LOSVD at each a is convolved to the thermal dispersion, and the results are integrated over the disk. 
No correction is made for a transition to e = at the center. 



Line profiles for various models are shown in Figure 3. The line shape appears at first glance very 
sensitive to the polytropic index n. But much of this sensitivity is a simple scaling; for the same central 
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Fig. 4. — A sequence of line profiles for the n = 3, e = 0.7 model, shown {left to right) as the disk precesses through a 
complete period. Gaussian noise has been added to simulate observations. The true mean velocity is zero, but continuum 
misidentification could lead to spurious measurements of velocity offsets. 



mass and limiting radius, disks of larger n show broader line profiles simply because they are more centrally 
concentrated. One also sees that the lines can be strongly asymmetric, even for eccentricities as small as 
0.1. Where the profiles are double peaked, the stronger peak arises from the side of the disk containing the 
apocenters of the orbits, and consequently can appear at either positive or negative velocity. For eccentricities 
of 0.5 or larger, the secondary peak can disappear entirely. Even though the line centroid remains at zero 
velocity relative to the central mass, in the presence of noise the asymmetric primary peak could easily be 
mistaken for a single narrow line with a large velocity offset. Figure 4 shows a sequence of line profiles 
through a complete precession period for the n = 3, e = 0.7 disk, with noise added. In principle, the slow 
time evolution of the line profile could be detectable for some types of objects. For a cold protostellar disk, 
equation (11) implies a precession period of 
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where the factor C = (4/n)(l — e^)^/^ is of order unity. For a galactic nuclear disk, 
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Offset or asymmetric emission lines are not uncommon in active galactic nuclei (AGN), but are usually 
qualitatively different from those shown above. Systematic offsets, e.g., of the [O III] line (Whittle 1985), 
are consistently toward the blue, and it is the line centroid, rather than the peak, that is shifted. Offsets of 
this sort should more naturally arise from outflows rather than orbital motion. AGN showing double-peaked 
lines characteristic of unresolved Keplerian motion may be modeled by axisymmetric relativistic disks when 
the blue peak is stronger than the red, the asymmetry arising from Doppler boosting (Chen, Halpern & 
Filippenko 1989). However, in a fraction of radio-loud AGN with disk-like Ha emission, the red peak is 
stronger than the blue. Eracleous et al. (1995) appeal to eccentric disks to model these objects. Their 
models, like those presented here, have constant eccentricity, but the line emissivity is taken to be a function 
of r rather than a. Eracleous ct al. (1995) consider general relativistic effects and estimate that differential 
relativistic precession would ruin the coherence of elliptic disks in AGN on tiinescales as short as decades, 
unless they were narrow rings. Bao et al. (1996) show that narrow rings can persist even in a strong field, 
and calculate the resulting line profiles. But it seems plausible that relativistic versions of the models of 
this paper could be found in which combined relativistic and hydrodynamic effects could drive a uniform 
precession, extending the lives of wide disks. 
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5. Discussion 

I have presented the properties of a particularly simple family of idealized models for eccentric Keplerian 
fluid disks, in which the internal pressure gradient drives a coherent precession of the apses of the constant- 
eccentricity streamlines. In the frame that rotates with the precession frequency, the pressure gradient 
balances the Coriolis force; the disks are therefore in global geostrophic balance, and the fluid motion is 
Keplerian, to first order in the pressure, in the rotating frame. 

In the limit of small eccentricity, the models can be viewed as perturbed axisymmetric disks, and their 
elliptic distortions interpreted as traveling acoustic waves. These waves are related to the "slow" m = 1 
modes examined by Trcmainc (2001), but arc not physically the same. The elliptic distortions have radial 
wavenumber fc = 0, which is maintained by a continuous refraction of the waves around the center, and 
therefore owe their existence to the particular density profile of the disk. The disk's self gravity, which 
creates the true standing waves studied by Tremaine, is ignored in the present case. Including finite disk 
gravity may restrict the slow acoustic waves to a limited range of radii and/or turn them into standing waves 
with k^O. 

One should keep in mind that these models are highly idealized, and neglect at least two important 
eff'ects present in real disks. First, fiuid elements in a three-dimensional disk will experience a varying vertical 
compression around every orbit, since the vertical {z) component of the Kepler force near the disk plane is 
proportional to z/r^. The fractional magnitude of the variation depends on the eccentricity and not on the 
disk thickness. Second, shear viscosity is known to be a controlling factor in true accretion disks, and is 
likely to affect the structure of larger-scale disks in galactic nuclei which may not be actively accreting. 

Turbulent viscosity can arise from any of a number of local instabilities, the leading candidate being the 
magnetorotational instability (Balbus & Hawley 1991) when magnetic fields are present. Angular momentum 
is transported efficiently by magnetic torques, and the instability's growth time is of order the orbital period. 
The instability will therefore affect the structure of any disk with even a weak magnetic field, as long as 
the ionization is suflaciently high that the flux is effectively frozen in the fluid. On the other hand, the 
instability may be suppressed in protostellar (Wardle 1999) or protoplanetary (Reyes-Ruiz 2001) disks of 
very low ionization. 

A purely hydrodynamical instability unique to non- axisymmetric disks has been studied by Goodman 
and collaborators (Goodman 1993; Ryu & Goodman 1994; Ryu, Goodman, & Vishniac 1996). The "eccentric 
instability" is intimately connected with the inertial oscillations, force-free epicyclic motions supported by 
any rotating or shearing medium. Because the spectrum of the inertial oscillations is continuous, there is 
always a mode of the right frequency available to be parametrically amplified by a periodic disturbance, such 
as a traveling acoustic wave. The instability is intrinsically three-dimensional — the motions of the amplified 
mode are inclined to the disk plane — and therefore would affect any finite-thickness analogue of the present 
models. 

Turbulence may also bo associated with vortices, whoso importance in disks is considered by Adams 
& Watkins (1995). They point out that, according to Kelvin's circulation theorem, vorticity should not 
be spontaneously created when pressure and density contours coincide and the viscosity is zero. Thus the 
models, as defined here, should not be unstable to vortex generation. Real disks, however, may be radiatively 
heated by their central objects and radiatively cooled by line or continuum emission. The temperature could 
then vary significantly around streamlines, in which case vortices could be generated. A careful treatment 
of the radiative transfer within the models would be needed to determine if this is the case. 
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Fig. 5. — Illustrating why a circular disk is not the minimum energy configuration for a given angular momentum. A pair of 
circular rings {filled circles) is plotted in the Lindblad diagram for a Kepler potential. Smooth curve shows the circular orbits. 
Open circles show the rings displaced so as to conserve total angular momentum and lower the total energy (IAE2I < 
leaving the outer ring eccentric. 

It may be tempting to argue that viscous dissipation will naturally circularize any disk on the grounds 
that, at a fixed angular momentum, the circular configuration is the state of lowest energy. However, 
this latter statement is true only for single orbits, and does not apply to disks where energy and angular 
momentum can be traded between orbits. To illustrate, consider the toy problem of two circular rings of 

equal mass. Figure 5 shows a Lindblad diagram for this system, with both rings (filled circles) lying on 
the locus of circtilar orbits (smooth curve) for the Kepler potential. Now imagine shifting the inner ring 
to a more tightly bound circular orbit (lower open circle), lowering its energy by AEi and its angular 
momentum by AL. Angular momentum can be conserved by depositing the same AL in the outer ring, 
along with some amount of energy Ai?2. Because of the shape of the circular-orbit locus, however, the shifted 
outer ring (upper open circle) can move away from the circular orbits (i.e., can become eccentric) even if 
\AE2\ < \AEi\. Thus, for any circular disk, there is an eccentric disk with the same angular momentum and 
lower energy. 

Energy arguments alone are therefore not sufficient to determine whether eccentric disks can survive in 
the presence of shear viscosity. A number of more detailed calculations have been attempted, but the results 
are still ambiguous. An analytic and numerical study of the viscous evolution of pressure-free eccentric gas 
disks and streams by Syer & Clarke (1992) finds that the eccentricity nearly always increases or remains 
constant under the influence of viscous effects. This result is corroborated by the purely anal}rtic work of 
Lyubarskij, Postnov, & Prokhorov (1994), who show that evolution at constant e is one of many possible paths 
for the disk. Ogilvie (2001) finds, in contrast, that with gas pressure included, eccentric motions dissipate 
on a viscous time scale. In three-dimensional smooth-particle hydrodynamic simulations (Mastrodemos & 
Morris 1998), broad, persistently eccentric disks are fotmd to form spontaneously in accreting winds. A more 
careful analysis of the simulations would be needed to clarify the detailed structure of the disks and the role 
of pressure in maintaining them. Disks whose radial profiles are determined by accretion processes may not 
evolve toward the special profiles required by the models in this paper, and consequently may be unable to 
synchronize their precession rates and remain eccentric. Btit the appearance of long-lived eccentric disks in 
simulations where no special effort is made to create them suggests that Nature may have ways to make 
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similar disks in abundance. 
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Appendix A. Surface Density of an Unaligned Keplerian Disk 

Equation (4) can be generalized to an unaligned planar disk in which the eccentricity and longitude of 
pericenter are given by functions e(a) and w{a) respectively. The result is 

i:(aE) = ^^ (^-^"y^" (481 

^' ' 27ra l-e2-a[e'(e + cos£;)+t37'e(l-e2)i/2sinS]' ^ ' 

where the eccentric anomaly E is measured from pericenter on each orbit. Alternatively, one might prefer 
to have the surface density in terms of a fixed polar coordinate system (r, (p). In this case, 

fija) (1-6^)1/2 

27ra 1 — — e'[2ae + r cos(0 — w)] — rew' sm{(f) — zu) ' 

where a{r, (j)) can be found efficiently by iteration using the prescription 

a,+i = (1 - k)aj + kr \ \ , ^^2 " (^0) 

Values of the relaxation factor k « 0.5 usually give rapid convergence. The above results generalize the 
formulae of Borderies, Goldreich, & Tremaine (1986) to arbitrary eccentricity, but remain valid only when 
the criterion for non-crossing of orbits (Ogilvie 2001), 

(e + ae'f + {aew'f < 1, (51) 

is satisfied. 



Appendix B. Helpful Integrals 

The following integrals are useful in evaluating the average precession rates in Section 2.3: 
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